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Abstract 

We show that a hybrid (nuclear + quark; matter) star can have a mass- 
radius refationship very similar to that predicted for a star made of purely 
nucleonic matter. We show this for a generic parameterization of the quark 
matter equation of state, and also for an MIT bag model, each including 
a phenomenological correction based on gluonic corrections to the equation 
of state. We obtain hybrid stars as heavy as 2 Mq for reasonable values 
of the bag model parameters. For nuclear matter, we use the equation of 
state calculated by Akmal, Pandharipande, and Ravenhall using many-body 
techniques. Both mixed and homogeneous phases of nuclear and quark matter 
are considered. 



1 Introduction 



It has long been hypothesized that some compact stars might actually be "hybrid 
stars" containing cores of quark matter. The observationally accessible features 
of compact stars include their mass and radius, and there have been various 
investigations of how the presence of a quark matter core would affect the mass- 
radius relationship of a compact star. The general conclusion has been that quark 
matter softens the equation of state, so that hybrid stars are predicted to have 
lower maximum mass (M < 1.7 Mo) than nuclear matter^ stars (M < 2.2 Mq) 
[Ill21l21llllilinilZllH]- On the observational side, all mass measurements are currently 
compatible, at the 2a level, with M < 1.7 Mq jSlE]- However, some are near the 
limit of compatibility, for example pulsar J0751+1807 in a white dwarf-neutron star 
binary system, whose mass is currently measured at 2.1(4) Mq ^T]. The error bars 
on these measurements will decrease over time so quark matter cores may seem to 
be on the point of being ruled out. Our purpose in this paper is to show that in fact 
it will be harder to rule out quark matter via M{R) observations than these simple 
considerations indicate. 

In an earlier paper [3] , we performed mass-radius calculations for hybrid and pure 
quark matter stars using a simple MIT bag model equation of state. In that model, 
correlations arising due to quark-quark interactions were neglected. However, the 
effects of the strange quark mass and corrections to the equation of state due to the 
pairing energy associated with color superconductivitv|12j . the formation of quark 
Cooper pairs, were incorporated. We found that color superconductivity boosts 
the pressure of the quark matter relative to nuclear matter, lowering the transition 
density (at fixed bag constant), but that the maximum mass was similar to that 
obtained in other work, namely about I.QMq. We found that in order to form stars 
near this upper bound it was necessary to set the bag constant to a low value, so that 
quark matter is very nearly stable, and the NM— ^QM phase transition occurs below 
nuclear saturation density, rigat = 0.16 fm~^: the heaviest stars consisted almost 
entirely of quark matter, with only a thin crust of nuclear matter. 

1.1 Overview of this study 

In this paper we will again use equations of state based on the MIT bag model for 
the quark matter, but we will include an additional parameter that imitates the 
effect of including perturbative QCD corrections. 

For the first stage of this analysis (Sec. E)) we will actually not use any model 
at all. We write down a purely phenomenological quark matter equation of state, 

^In this paper we only consider two-flavor nuclear matter. Introducing hyperons or kaon 
condensation into the nuclear matter also softens the equation of state and lowers the maximum 
mass. p. 
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consisting simply of a power series expansion in the quark chemical potential, /i 



where 04, 02 and B^s are independent of jj. We will show that for 04 ~ 0.7 
(which we will see later is physically reasonable), one can obtain heavy hybrid 
stars (M ^ 2 Mq) while still ensuring that the NM— >QM phase transition occurs 
above nuclear saturation density. In fact, we show that it is possible to mimic the 
mass-radius behavior of nucleonic stars over a wide range of masses. 

We then go on to the second stage of our analysis, in which we use a quark 
matter equation of state based on a physical model: competition between a "normal" 
unpaired quark matter phase and the CFL color-superconducting phase, in an MIT- 
bag-model formalism. This corresponds to giving 02 and Bes a simple step-function 
dependence on the chemical potential, with the step occurring at the transition 
between these phases. The details are determined by the microscopic parameters of 
the model, the strange quark mass and pairing gap A. 

The coefficient 04 = 1 — c is a rough parameterization of QCD corrections to the 
pressure of the free-quark Fermi sea, and previous calculations ^21 show that the 
value 04 ~ 0.7, c ~ 0.3 is reasonable. The physical model also allows us to calculate 
the behavior of charged phases, and hence to study inhomogeneous "mixed" phases 
of nuclear matter and quark matter. 

The results from the quark matter model turn out to be very similar to those 
obtained with the simple parameterization. We also find that, given our current 
(and, likely, future) ignorance of the high-density values of basic parameters like 
the strange quark mass and bag constant, there are no characteristic features 
of the M{R) relationship that could be used to verify the presence of color 
superconductivity in the quark matter core. 

Our conclusion is that the maximum mass of hybrid stars is mainly determined 
by the size of the QCD corrections to the coefficient 04 in the quark matter equation 
of state (jl.ip . and that for reasonable values of hybrid stars can be as heavy 
as 2 Mq. This result is robust: it is not affected when we move up to a more 
sophisticated model, nor by the introduction of mixed phases. There is therefore 
little reason to expect it to change when the model is made even more complicated, 
e.g. by including kaon-condensed [H], crystalhne [12], mixed [12] or gapless [T7| ITH] 
phases, or allowing continuous //-dependence in the strange quark mass or pairing 
gap (see below). 

1.2 Other approaches 

Our calculations use a basic MIT bag model, in which the bag constant and quark 
masses are assumed to be density-independent. Other approaches are certainly 
possible. Within the MIT bag model one can use a density-dependent bag constant, 
although this does not appreciably change the maximum mass prediction [0]. The 





02 + BcS, 
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density-dependence of the constituent quark masses and the color superconducting 
gap A can be estimated by using an NJL model instead of a bag model, with coupled 
mean-field Schwinger-Dyson equations for the masses and gaps. Such models give a 
high effective bag constant and quark masses, and typically predict small numbers 
of strange quarks, with 2SC rather than CFL color superconductivity, but again the 
maximum masses turn out to be of order 1.6 Mq [H [7], although folding certain 
Gaussian form factors into the 4-fermion interaction can give masses up to 1.8 Mq 
[T^ . While these models are well motivated theoretically, their specific predictions 
relating to the density dependencies of quark masses and the effective bag constant 
remain untested. In this work we adopt a minimal approach and retain quartic and 
quadratic powers of the chemical potential in our expression for the free energy. This 
will allow us to do a parameter study independent of any specific model. It would 
be interesting to see whether including a reasonable estimate of QCD corrections in 
the NJL model increases the maximum mass in that context also. 

2 The nuclear matter equation of state 

Our treatment of nuclear matter is completely standard: at densities above half 
nuclear saturation density (rigat) we employ the equation of state of Akmal, 
Pandharipande, and Ravenhall (APR) [201 • At lower densities we use the standard 
tabulated low-density equation of state j2Ill22|- Our previous studies of hybrid stars 
12S] used the relativistic mean field Walecka model [25 to describe nuclear matter 
or APR for beta stable charge-neutral nuclear matter. The relativistic mean field 
is an effective description of the nuclear matter that is constrained by properties of 
nuclear matter at saturation density. The mean field approximation ignores many- 
body correlations that could play an important role. APR use the variational chain 
summation (VCS) to include these correlations in calculating the equation of state 
of nucleon matter fl^ 1^ . They employ a realistic non-relativistic Hamiltonian 
with the Argonne fig [27j two-body potential and the Urbana IX [2H] three nucleon 
interaction. An equation of state as a function of baryon density and proton fraction, 
Xp is obtained by interpolating between the pure neutron matter (PNM) Xp = 
and symmetric nuclear matter (SNM) Xp = 0.5 results using a generalized Skyrme 
interaction containing momentum and density dependent delta function interactions 
[2ni described below. 

The energy is evaluated for a variational wave function which takes into 
account many-body correlation effects. It is composed as a symmetrized product 
of two-body correlation operators, Fij acting on the Fermi gas wave function. 
The correlation operators are written as a sum of terms which include operators 
appearing in the Hamiltonian. The two-body cluster contribution to the energy is 
minimized by Fij which satisfy the Euler-Lagrange equations determined within this 
ansatz. Heuristically, two-body operators which appear in the Hamiltonian induce 
correlations between particles whose spatial dependence is approximated by solving 
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a two-body Schrodinger-like equation subject to suitable boundary conditions. 
The variation of the wave function is effected at the two-body level by varying 
parameters appearing in this equation. Though this wave function neglects three- 
body correlations it is estimated to be accurate to a few MeV/nucleon in SNM 
and about one MeV/nucleon in PNM at nuclear density. This accuracy is achieved 
through the inclusion of many-body effects via the VCS technique. Comparable 
accuracy is obtained with the other many-body techniques like the Brueckner-Bethe- 
Goldstone method of Ref . PU] . 

The APR equation of state exhibits a transition from a low density phase (LDP) 
to a high density phase (HDP) having spin-isospin order, possibly due to neutral 
pion condensation, in PNM at a density of ~ 0.20 fm~^ and in SNM at ~ 0.32 fm~^. 
The VCS calculations of the energy of PNM and SNM are extrapolated to general 
values of Xp using a function of the form 

eiv(P, Xp) = h /(p, Xp) \ Tp + h /(p, 1 - Xp) \ Tn 



2m/vr / \2mjv 

+ g{p, Xp = 0.5)(1 - (1 - 2xpf) + g{p, Xp = 0)(1 - 2xpf (2.1) 

motivated by a generalized Skyrme interaction. Here e^v is the total nuclear energy 
density, r„.p are the neutron and proton Fermi gas kinetic densities, and f{p,Xp) 
and /(p, 1 — Xp) are functions which parameterize the effective mass of the nucleons 
and g{p,Xp = 0) and g{p,Xp = 0.5) are potential energy terms. These functions 
are parameterized to fit the energies of PNM and SNM calculated in VCS. Separate 
parameterizations are used in the LDP and HDP for the functions g. The energy 
density of nuclear matter will be used to determine the allowed equation of state of 
hybrid stars for the case of a sharp transition to quark matter and to determine the 
allowed phases for a mixed transition to quark matter. 



3 A simple phenomenological quark matter equation of 
state 

One of the main points of this paper is to show the effects of including a parameter 
in the quark matter equation of state that roughly incorporates the effects of gluon- 
mediated QCD interactions between the quarks in the Fermi sea. We first do this in 
the context of a simple parameterization of the quark matter equation of state, and 
show that the resultant hybrid stars can have mass-radius relations very similar to 
that of pure nuclear stars, with mass up to 2 Mq. Later (Sec. |3)) we will show that 
these conclusions remain true in a more sophisticated model. 

The simple phenomenological parameterization of the quark matter equation of 
state is 

3 3 

^QM = 04 + p' + 5efr , 

a4 = 1 — c , 
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where the parameters 04, 02, Bcs are independent of yU.. 

1) The quartic coefficient 04 = 1 — c. 

For quark matter consisting of 3 flavors of free non-interacting quarks, c = 0, so 
a4 = 1 (see the discussion of the physical model of quark matter in SecE]). However, 
once QCD corrections are taken into account, we expect c 7^ 0. The QCD corrections 
to the quark matter equation of state were ffist evaluated to 0(«s) by Freedman 
and McLerran [3T]; then Fraga, Pisarski and Schaffner-Bielich ^2] (FPS) showed 
that the 0{al) pressure for three massless flavors can be approximated by a bag- 
model-inspired form, 

^-i ~ ^ (1 - /i' - 5eff . (3.2) 

Matching to the 0(«s) perturbative calculations in the interval /i ~ 300 — 600 MeV, 
they found that B^s varies widely with renormalization scale (see also Ref. |32j). 
but they consistently find c ~ 0.37. We do not use FPS's specific values of B^s and 
c because, as they observe, the QCD coupling is strong at the density of interest for 
compact star physics, so there is no reason to expect the leading order calculation 
to be accurate. However, we take their results as indicating that QCD corrections 
are not negligible, so one should include c as an additional parameter in the quark 
matter equation of state, with a value of order 0.3. 

2) The quadratic coefficient 02. 

As we will see in Sec. lU the term can arise from the strange quark mass (which 
increases the free energy) or color superconductivity (which reduces it). If chiral 
symmetry remains broken in the light quark sector due to a robust (qq) condensate, 
then the large (~ 100 — 300 MeV) constituent quark masses of the up and down 
quarks would also result in fi"^ term similar to that due to the strange quark mass. 
For color-flavor-locked quark matter, 02 = m"^ — 4A^. For now, we simply include 
a2 as a phenomenological parameter. 

3) The bag constant and the transition density pc- 

In our parameterization the effective bag constant simply accounts for the free energy 
contribution that is independent of fi. While this is related to the vacuum pressure 
its numerical value in our parameterization need not be the same as in early bag- 
model studies of hadron phenomenology . The effective bag constant is unknown 
and difficult to calculate or measure, so we will not use it as a parameter, since that 
would obscure the fact that part of the effect of varying other parameters may simply 
be a renormalization of the unknown parameter Bcr. To expose the physically 
significant effects of varying 04, 02, Beg we will specify a more physical quantity, 
the maximum density pc of nuclear matter, i.e. the density at which the NM— >QM 
transition occurs. The structure of the star is then calculated as a function of 04, 02 
and Pc- 

Two subtle points arise in such a reparameterization. Firstly, as we will see 
below, for values of the perturbative correction parameter c around 0.3, the quark 
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matter and APR nuclear matter equations of state have almost exactly the same 
shape over a wide range of pressures. This can lead to multiple phase transitions 
back and forth between NM and QM. Of course, when the two phases have almost 
identical equations of state, it does not matter (for mass and radius calculations) 
where transitions between them occur. We therefore simply choose pc to be 
associated with the location of the first transition. 

Secondly, when we use the more complicated model in Sec. 01 we will allow 
for the possibility of mixed phases, which blur out the NM— i>QM transition over a 
range of densities and pressures, making it hard to identify "the" transition density. 
However, for the purpose of fixing the bag constant we do not have to allow mixed 
phases. We will therefore define B^sia^, a2, pc) as the value of the bag constant that 
would give a sharp NM— >-QM transition at a nuclear matter density if only charge 
neutral bulk phases were permitted (as would happen if the NM-QM surface tension 
were ^ 40 MeV/fm^). 



3.1 The physical effects of the "perturbative correction" parameter 

We now discuss the physical importance of the QCD correction c. Firstly, it is 
clear analytically that, at fixed bag constant, c has very little effect on the relation 
E{p) between energy density and pressure for quark matter, which enters into the 
Tolman Oppenheimer Volkoff (TOV) equation [33]. This can be seen by setting 
02 = (i.e. neglecting quark masses and pairing), in which case 

P= (l-c)4^-5eff 

=^ E = 3p + 4i?cfr independent of c 

However, this does not mean that c is unimportant. Clearly c makes a dramatic 
difference to p{p), so it strongly affects the position of the NM— >-QM transition. 
Moreover, as described above, we are not working at fixed Bcs: when we change c 
we keep the transition density pc fixed, with a resultant change in Beg. 

We illustrate the effect of non-zero c in Fig. ^ which shows p{p) and E{p) for 
APR nuclear matter and for the phenomenological description of quark matter, with 
a2 = (150 MeV)^, for the cases c = and 0.3 (tuning B^s to keep the transition 
density p = l.Srisat)- We see that the c = 0.3 quark matter equation of state is 
very similar to APR over the pressure range 10 to 200 MeV/fm'^. In fact, close 
examination of the p{p) relationship shows that there are three phase transitions 
with increasing density: from APR nuclear matter at low density to quark matter, 
then back to APR, then back to quark matter again. On this basis, we expect stars 
containing quark matter with c ~ 0.3 to show M{R) curves very similar to those of 
nuclear matter stars, making them correspondingly difficult to rule out from M{R) 
measurements alone. In Sec. 13.31 we will see that this is indeed the case. 
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)i (MeV) Pressure (MeV/fm^) 



Figure 1: Equations of state for APR nuclear matter (solid line) and for quark matter 
with the phenomenologically parameterized equation of state (|;-{.lj) with 02 = (150 MeV)^ 
(dashed curves: upper curve has no perturbative correction (c = 0), lower one has c = 0.3). 
For each quark matter equation of state the bag constant B^s was fixed by requiring that 
nuclear matter give way to quark matter at pc = l.Sngat- The figure shows that when we 
include perturbative-type corrections, quark matter has a p{p) relation almost identical to 
that of APR nuclear matter, and its E(p) is also very similar over the range of pressures 
that is relevant to compact star masses. 



The p{fi) relations for c = 0.3 quark matter and APR are so similar that it can 
be difficult to say where the phase transition really occurs. In Fig. |21 we show 
the pressure difference for the parameters of Fig. ^ {bH^'^ = 141.2 MeV, solid 
curve) and for a higher value of the bag constant (-B(!g^ = 141.5 MeV, dashed 
curve). At -Bgg^ = 141.2 MeV, the three transitions are clearly visible. We see 
that when the bag constant is raised to -Bg^^ = 141.5 MeV to obtain a transition at 
higher density, the value of the chemical potential at the transition, and hence the 
density, jumps discontinuously from /i = 338 MeV (p = l.Qrisat) to /i = 458 MeV 
(p = 4.3nsat)- This is illustrated in the right panel of Fig. |21 It is therefore 
technically impossible to choose Bcs so as to obtain a (first) transition density 
in the range 1.9nsat < Pc < 4.3rasat- 

It is very interesting to speculate on the possibility of multiple phase transitions 
inside compact stars, but this feature is highly sensitive to the precise relative shapes 
of our quark matter and APR equations of state. We cannot claim to know these 
to the level of accuracy (a few percent, i.e. a few MeV/fm^) that would be required 
to say whether the number of transitions is 1, 3, or even 5. Our main message 
in this paper is that for reasonable values of the quark matter parameters the 
equations of state may be very similar. A realistic interpretation of Fig. |21 is that 
at -Bgg^ = 141.5 MeV the critical density jumps rapidly from l.Qrisat to 4.3nsat- 
Any Pc in that range therefore corresponds to -Bg^^ = 141.5 MeV. (Pictorially, this 
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Figure 2: The pressure difference between quark matter and APR nuclear matter for 
04 = 0.7, 02 = (150 MeV)^, illustrating how the near-identity of their relations leads 
to difficulties in defining a unique transition density. In the left panel, the solid curve is 
for -Bgg"^ = 141.2 MeV and shows three transitions, the first occurring at ^ = 324.4 MeV, 
corresponding to pc = l.Snsat (Fig-©- Raising B^s gives a transition at higher ^, but it 
is technically impossible to achieve a (first) transition in the range 338 < fi < 458 MeV 
(1.9nsat < P < 4.3nsat): at -Bg^^ = 141.5 MeV (dashed curve) the transition jumps from 
the bottom of this range to the top. The right panel illustrates how this translates to a 
discontinuous behavior in dependence of pc on B^s- 

corresponds to blurring out the curves in Fig. m by a few MeV/fm^) This is the 
procedure we will follow in determining Sgfr for given p^, mg, A when we investigate 
how the maximum compact star mass depends on these parameters. 

One can extrapolate from Figs. Q and |21 to predict what will happen for even 
larger values of c. For suitable values of the bag constant there will be a small range 
of densities (near nuclear density) where quark matter has higher pressure than 
APR nuclear matter, then at intermediate densities APR will be favored again, and 
finally quark matter will win at high densities. When such equations of state are 
used to construct compact stars, the result will be a star containing a shell of quark 
matter with APR nuclear matter outside and inside it. This also leads to large 
masses, but only because most of the star is nuclear matter. In the rest of this 
paper we will discuss values of c up to 0.3, which is suggested by Fraga et. al.'s fit 
to the two-loop equation of state ^2] , and also offers the possibility of heavy stars 
with a considerable quark matter fraction. 

3.2 Mass-radius relationship: a hybrid star that "looks nuclear" 

Before surveying a wide range of values of the parameters 04 = 1 — c, 02, Pc, we first 
display the M{R) curve for some specific cases, showing how masses approaching 
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APR + Phenomenological QM EoS 




R (km) 

Figure 3: M{R) relationship for hybrid stars involving quark matter obeying the 
phenomenological equation of state p.lf) with 02 = (180 MeV)^. In each case the bag 
constant was chosen to give the desired transition density. The line with black dots is the 
M{R) relation for a pure nuclear APR star. All the phases are neutral and homogeneous. 
We see that the maximum mass Mmax is very sensitive to the QCD correction c, but not 
to the transition density. For c = 0.3, which is close to the value suggested in Ref. |13j . 
-^max ~ 1-9 Mq. Note also that the c = 0.3 equation of state is so similar to APR that 
its M{R) curve is almost identical to that of pure APR. 

2 Mq can be achieved. The M{R) relation is obtained by solving the TOV equation 
in the standard way, as described in Ref. [3]. We choose a range of central pressures, 
and integrate the TOV equation outward until the pressure drops to zero, which 
marks the surface of the star, and the integrated energy density yields the mass jHS] • 
In Fig.iniwe show the resulting M{R) curves for quark matter with a2 = (180 MeV)^. 
We performed calculations with c = and with c = 0.3 ^3]- We tuned the bag 
constant to give a range of homogeneous neutral-matter transition densities Pc = 2 

to 4nsaf 

The most obvious feature is the dependence on c: compared to the stars with 
c = 0, hybrid stars containing quark matter with larger perturbative correction 
c = 0.3 are significantly heavier and a little larger, with masses approaching 2 Mq, 
and radii around 11 km. It is also striking that stars with c = 0.3 have an M{R) 
relation essentially identical to that of APR nuclear matter, up to masses around 
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APR + Phenomenological QM EoS 




R (km) 

Figure 4: The radius at which the first (lowest-density) transition from nuclear matter to 
quark matter occurs, for the families of stars whose M{R) relations were given in Fig. 13 
Note that this radius is quite large for the heavier stars. In the c = 0.3 case there may be 
transitions back and forth between quark and nuclear matter at higher density (smaller 
radii) . 

1.9 Mq. This is because quark matter with c = 0.3 has an equation of state very 
close to that of APR nuclear matter over the relevant range of pressures. This was 
discussed in Sec. 13.11 

In Fig. m we show, for the same family of stars, the radius at which the first 
(lowest-density) transition from nuclear matter to quark matter occurs. Note that 
this radius is quite large for the heavier stars. For values of c close to 0.3 there are 
transitions back and forth between quark and nuclear matter as one goes deeper 
into the star, because the equations of state are so similar. While such multiple 
transitions cannot be ruled out a priori, we suspect that they are not physical. They 
are an artifact of our model description of the nuclear and quark phases and not a 
robust prediction. Here and in the rest of the paper we will ignore this possibility 
and entertain only one transition from nuclear to quark matter, beyond which we 
will use the quark matter equation of state. While this will make little difference to 
the structure of the star, since the nuclear and quark matter equations of state are 
both very similar in this regime, it would have important consequences for transport 
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APR + Phenomenological QM EoS 




sqrt(a2) (MeV) 

Figure 5: Dependence of Mmax on 02 for hybrid stars involving quark matter obeying the 
phenomenological equation of state Increasing 02 decreases Sefr at fixed pc, giving 

rise to larger maximum masses. We show results for two transition densities, pc = 3nsat 
and /9c = 2 ngat • The maximum mass is not very sensitive to the transition density. 

properties. 

3.3 Maximum mass as a function of quark matter parameters 

Having seen that with a perturbative correction c set to a reasonable value we can 
increase the maximum hybrid star mass, we now look at how maximum hybrid star 
mass depends on the parameters of our phenomenological equation of state ()3.1|) : 
04, 02, and the transition density pc- 

In Fig. El we plot the maximum mass, Mmax obtained by varying the central 
pressure, and choosing the heaviest resulting stable star that contains some quark 
matter. We calculated Mmax as a function of 02, for 04 = 1, 0.8, 0.7, 0.65. In each 
case we repeated the calculation for two different transition densities, pc = SrZsat 
(solid lines) and pc = 2nsat (dashed lines). The lines in Fig. El end when there is 
no longer a stable hybrid star. This corresponds to the hybrid star branch in the 
M{R) plot having a positive slope along its whole length. 

We see that masses up to 2 Mq can be obtained by choosing a small value of 04, 
and setting ^/a2 to an appropriate value in the range 150 MeV to 250 MeV. (Recall 
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that in unpaired quark matter ^Ja^ corresponds to the strange quark mass.) This 
is understandable since the effective bag constant must decrease with increasing 
02 at fixed p^. it is well known that in stars dominated by quark matter the 

M^ax ~ 1/v^ EH- 

One might want to ask how large the quark matter core is for the heaviest 
stars, but this question does not have an easy answer. For 04 ~ 0.7 (c ~ 0.3), the 
quark matter and nuclear matter equations of state are so similar that the transition 
density is not precisely defined (as we saw in Fig. 12)) so it is not clear where the 
transition occurs: there is a whole range of transition densities that correspond to 
the same bag constant. This is why, for lower values of 04, the maximum mass is 
not very sensitive to the transition density. Fig. |3 shows that the first transition 
from nuclear to quark matter occurs at quite large radii for the heavy stars. 

3.4 Conclusion 

We have seen that the very simple phenomenological parameterization p.ip of the 
quark matter equation of state allows us to survey the main effects of the various 
parameters, and makes it clear that 04 = 1 — c is the most important in determining 
the maximum mass. With a value corresponding to free quark matter (c = 0) we 
find a maximum hybrid star mass of about 1.6 M©, similar to the value found in 
previous studies. With a value corresponding to QCD corrections of a plausible 
strength (c ~ 0.35) we find a maximum mass of about 2 Mq. 

In the next section we show that this conclusion is robust. We repeat our analysis 
with a more complicated quark matter equation of state, based on the expected 
physics of dense quark matter. We allow competition between paired and unpaired 
phases and mixed quark and nuclear matter phases. We find that these factors have 
no effect on our essential findings regarding the maximum mass of hybrid stars. 

4 A more sophisticated approach: a physical model of quark 
matter, and mixed phases 

The model we will use in this section is more complicated and physically well- 
founded than the phenomenological parameterization of Sec. El but it is still a 
relatively simple MIT bag model. We will allow competition between two phases of 
quark matter, the three-flavor unpaired phase and the color-fiavor-locked (CFL) 
color-superconducting phase. The phase with the lowest free energy (highest 
pressure) is favored. 

The parameters of this model are the strange quark mass m^, the CFL pairing 
gap A, the density pc of the nuclear to quark matter transition (which determines 
the bag constant), and a QCD-inspired correction parameter c. This is similar to 
the phenomenological parameterization of Sec. IHl but the effective parameter 02 is 
now calculated within a model. 
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At asymptotically high chemical potential on the order of 10*^ MeV, the CFL 
phase is known to be the ground state of 3-flavor quark matter, though it may 
persist to much lower chemical potential. At lower densities the strange quark mass 
becomes non-negligible relative to the chemical potential, and at ml/ ^ = 4A there 
is a transition to unpaired quark matter ^7\. ^ When we include the nuclear matter, 
there is a three-way competition between nuclear matter, unpaired quark matter, 
and CFL quark matter. 

Now that we have a physical model in hand, we can calculate the pressure of 
charged as well as neutral phases of quark matter. This means that we can study 
mixed as well as homogeneous phases. The transition from nuclear matter to quark 
matter can proceed via a mixed phase [2H] , in which there is charge separation, and 
a positively charged nuclear phase interpenetrates with a negatively charged quark 
matter phase, yielding a globally neutral inhomogeneous phase. This will only occur 
if the surface tension at the boundary of the two phases is low enough. Otherwise 
there will be a sharp interface between the two homogeneous neutral phases. The 
critical surface tension is o"c ~ 40 MeV/fm^ |39| 1^ . Since the surface tension 
is completely unknown, we will separately consider both mixed phases and sharp 
interfaces in our calculations. 



4.1 Unpaired 3-fiavor quark matter 

In the naive bag model, where perturbative corrections are ignored, the free energy 
of free quark matter consists of the kinetic contribution from a degenerate free gas 
of three colors of relativistic quarks and the negative vacuum pressure from the bag 
constant, B^s 



^nnpil^u, IJ^d, l^s) = ^ J dp p^{^J 'p^ + mf - ^i) + BcS (4.1) 

i=u,d,s ^ 

where iiu = l^ — |/^e, /^d = /^s = /^ + |/^e are the individual quark chemical potentials 
and /i and /ie are the baryon and electron chemical potentials respectively. We 
may neglect the quark masses of the up and down quarks in quark matter since 
'flight ~ ~ ^ yU. On the other hand, the strange quark mass is not negligible 
compared to In this work we will assume that m<, < /i so that an expansion in 
powers of m^/ yU is meaningful and study three flavor quark matter. 

In neutral unpaired quark matter, the electron chemical potential is determined 
by the condition of charge neutrality and is given by 

2 4 / fi 

mt m: „ / m 



^Actually, before that, at ml = 2/iA, there is a transition to a gapless phase jT^, and at 
TOs ~ "^light^^^'^ there is the possibihty of A'*^ condensation in the CFL phase ^Jj: we ignore 
these additional comphcations because their contributions to the pressure are of order to^, which 
is just a renormahzation of the bag constant, which for us is a non-physical parameter. 
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Substituting this in Eq. 14. II and expanding in powers of rris/fi we obtain 



where we have dropped terms of order mf///^ and higher, and we have introduced 
a parameter c corresponding to the QCD- inspired corrections of Ref. just as 
we did in Sec. 0] From Eq. ()4.3p we see that the expansion in powers of m^/ fi is 
rapidly convergent, even for ~ /i. The contribution to the free energy from 
the m^/fi^ term is less than one part in 10^ of the /x'^ term when ~ 300 MeV 
and /i = 350 MeV. This means that neutral unpaired quark matter is really just a 
particular case of the phenomenological parameterization explored in Sec. El The 
EPS parameterization relies on the analysis of Ereedman and McLerran which is 
rigorous for massless quarks ^T]. We therefore apply the effects of the "perturbative" 
QCD correction only to the (i.e. massless) part of the free energy. Again, this is 
similar to the phenomenological parameterization, where the QCD correction is a 
modification of the coefficient of the /i^ term. 

If we do not impose the neutrality condition, but expand in powers of fie/ /i as 
well as mg/fi, we find that the expansion does not converge nearly as well, so for 
charged unpaired quark matter we must use the full form of the free energy Eq. (j4.1|l . 



4.2 Color-fiavor-locked (CFL) quark matter 

In color-flavor locked matter, the pairing locks the Eermi momenta of the all the 
quarks to a single value, requiring the number densities of up, down and strange 
quarks to be equal [lOj. This costs free energy, which is offset by the pairing 
contribution 

nA = -^ AV + (^(A')- (4.4) 

Calculations of A with effective interactions yield values in the range 10 — 100 MeV 
[T^ for /i in the range 300 — 600 MeV. So it is reasonable to retain only the leading 
order (in powers of A) contribution. As in the case of the strange quark mass, the 
A^ contribution has a weak (logarithmic) dependence on /i and its contribution to 
the equation of state is indistinguishable from B^s- 

When color and electric neutrality is imposed, there are no electrons in the CEL 
phase since there are equal numbers of up, down, and strange quarks, so the electron 
chemical potential /ig = 0. Expanding the free energy in powers of m^/ /i, we find[ 

neutral /i \,.4 , 'J''''sP p 



^CFl' quarks " ^^2^^ 4^2 ^2 

+ 
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^neutr. Sm^ - 48A V ( m[ 



where we have assumed that the unpaired and CFL phases have the same bag 
constant. 

We see from the first hne of this expression that the equation of state of neutral 
CFL matter is controlled by the strange quark mass and pairing gap in a combination 
that corresponds to the parameter 02 in the phenomenological parameterization 
fl3.H) , if we identify 

a2 = m1- . (4.6) 

Note that the terms do not follow such an identification, but they are irrelevant 
because they are /^-independent, so they are a renormalization of the bag constant, 
which we choose to give some specified transition density. Therefore the neutral CFL 
equation of state, to the order that we have expanded it, is rigorously a function 
of the QCD-inspired correction c, the chosen transition density pc, and — 4A^. 
However, the equation of state for neutral quark matter in general is not rigorously 
a function of — 4A^ only, because for a given and A there is competition 
between the neutral CFL phase and the neutral unpaired phase, whose equation of 
state is affected by rris alone. However, we expect the dependence on the linearly 
independent variable Ami + A^ to be very weak. By independently varying and 
A we have verified that this is the case. Of course, mixed phases are complicated 
by the fact that they involve charged quark matter whose equation of state depends 
on He as an additional parameter. 

We want to construct mixed phases so we must also know the equation of state 
for charged CFL matter with /ie 7^ 0. This has contributions from the quarks, 
the electrons, and from Goldstone bosons. The Goldstone bosons arise due to the 
spontaneous breaking of chiral symmetry, analogous to the mechanism in vacuum 
|41j . The quark contribution is independent of fie'- CFL-paired quarks form an 
insulator with gap A, so as long as /ie < A there are no charged quasi-quark 
excitations, so Q = dQ/dfj,e = 0. 

The Goldstone bosons are a consequence of spontaneous breaking of chiral 
symmetry jU |2S1 and like the octet of pseudoscalar mesons of QCD in vacuum, 
can be described by an effective chiral field theory jl21- When the electron chemical 
potential exceeds the mass of the lightest negatively charged meson, which in the 
CFL phase is the tt", these mesons condense ^3]. The free energy contribution of 
the meson condensate is 

where the parameters are |32] 

. (21-81n2)//2 2 3A2 ^ ^ ^ 

" 36^^2 ' = ^"^«("^« + ^^d) ■ (4-8) 
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APR + bag model QM EoS 




8 10 12 14 

R (km) 



Figure 6: M{R) relationship for hybrid stars involving quark matter obeying the simple 
bag model equation of state, with various values of c and p^, and = 180 MeV and 
A = 0. The line with black dots is the M[R) relation for a pure nuclear APR star. The 
solid lines are for hybrid stars with homogeneous neutral APR and CFL phases, the dashed 
lines are for hybrid stars with mixed APR+CFL phases. We see that the maximum mass 
-^max is very sensitive to the perturbative correction c, but not to the transition density, 
or the occurrence of mixed phases. For c = 0.3, which is close to the value suggested in 
[IHI, AW ~ 1.9 Mq. 

We used m„ = 3.75 MeV, = 7.5 MeV. Finally, the free energy contribution from 
electrons and muons is given by 

^leptons^^^) ^ ^ - y dpp\^p^ + mi-^^,). (4.9) 

i=e~ ^ 

The total free energy for CFL quark matter is then 

f^e) = ^CFL-UM + fi8FL(/i, /ie) + n'^^'^if^e) (4.10) 

4.3 A bag-model hybrid star that masquerades as a neutron star. 

Using the bag model equations of state for unpaired and CFL quark matter, we 
can proceed as in Sec. El to solve the TOV equation and obtain M{R) curves for 
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hybrid stars. We first show a specific example, in Fig. IHl We have chosen the 
bag model parameters to give the same equation of state for neutral matter as was 
studied using the phenomenological parameterization in Fig. El Now, however, we 
can also study mixed phases. The solid lines, for homogeneous neutral quark and 
nuclear matter phases, are the same as in Fig. 01 The dashed lines show the M{R) 
relation when mixed phases are allowed. We see that although the overall shape 
of the M{R) curve is quite different when mixed phases are present the maximum 
mass is not significantly affected. The maximum mass of a compact star is therefore 
insensitive to the surface tension of the interface between quark matter and nuclear 
matter. This is because the maximum mass configuration is characterized by a 
baryon density that is large compared to the transition density and most of the 
star is either in the homogeneous quark matter phase, or in a mixed phase that is 
dominated by quark matter. 

We can also see that, as in the phenomenological model explored in Sec. O the 
maximum mass is very sensitive to the QCD correction c, and relatively insensitive 
to the transition density. All the curves in Fig. (HI are for = 180 MeV, A = 0, so 
from this figure we cannot judge the sensitivity to those parameters. 

4.4 Maximum mass as a function of bag model parameters 

In Fig.[7|we give a more complete picture, by plotting the maximum mass of the star 
as a function of ^y^m'^ — AA'^. Varying the independent variable, ^^Aml + A^, has 
negligible effect, as one would expect. We plot this relationship for = 3nsat and 
Pc = 2r2sat- The stars with small — 4A^ correspond to hybrid stars containing 
CFL quark matter while those with large m^— 4A^ correspond to stars with unpaired 
quark matter. We see that the largest masses are obtained by turning — 4A^ up 
to a large value, corresponding to a strange quark mass in the range 200 MeV to 
300 MeV. This is the same situation that was described in Sec. l3.3l since the effective 
bag constant must decrease with increasing — 4A^ at fixed pc- As mentioned 
earlier, stars dominated by quark matter have Mmax ~ 1/ \/ij^ff pHj. We also see 
that the value of the transition density is not very important: Mmax is only weakly 
dependent on the transition density, particularly when the perturbative correction 
is included. 

5 Conclusions 

We have studied the mass-radius relationship for hybrid compact stars, with a 
nuclear matter crust (described by the APR equation of state) and a quark matter 
core. We used a phenomenological parameterization ()1.1|) of the neutral quark 
matter equation of state, and also a simple MIT bag model that allowed us to 
construct mixed phases. In both cases we included a QCD correction parameter 
c in the quark matter equation of state, and we found that increasing its value 
from zero (no QCD corrections) to reasonable value (c = 0.35 [T3|) increases the 
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APR + bag model QM EoS 




Figure 7: Maximum mass Mmax for hybrid stars containing quark matter obeying the 
simple bag model equation of state. We show the dependence on y^m^ — 4A'^. Increasing 
this quantity decreases i?efr at fixed pc, giving rise to larger maximum masses. The 
maximum mass is only weakly dependent on the transition density, particularly when the 
perturbative correction is included. 



maximum hybrid star mass from about 1.6 Mq to about 2 Mq. This is clear from 
Fig. El (for the phenomenological parameterization) and Fig. [7| (for the bag model). 
The reason is that increasing c hardens the quark matter equation of state, making 
it almost indistinguishable from the APR nuclear equation of state (see Fig. Q). It is 
important to note that we achieve these masses with reasonable transition densities 
Pc of order 2 to 3 times nsat- Our stars have a proper crust of nuclear matter: we are 
not resorting to low transition densities that yield "hybrid" stars that are actually 
quark stars with a tiny shell of nuclear matter around the outside. 

We conclude that it is harder than previously thought for a simple mass 
measurement to rule out the presence of quark matter in compact stars. Currently 
published measurements of the masses of compact stars are all consistent with 
a maximum mass of Mmax ~ 1-7 Mq at the 2a level. Our results show that 
maximum masses of up to 2 Mq can be accommodated by models of hybrid stars 
with reasonable quark matter equations of state. 

We also note that our M{R) curves are consistent with the constraint obtained 
from measurements of red shifts of Iron absorption lines in the low-mass X-ray 
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binary EXO0748-676 |13]. If that constraint were plotted in our Fig. IHl it would 
intersect our c = 0.3 curve at M ^ 1.7 Mq. 

This naturally raises two related questions: what sort of mass or radius 
observation would provide evidence against the presence of quark matter in neutron 
stars? And what sort of observation would provide evidence for the presence of 
color-superconducting quark matter in particular? 

Obviously a mass measurement above 2Mq would give reason to doubt the 
presence of quark matter. Other than that, it seems difficult to diagnose the presence 
of quark matter via M{R) measurements. The regions of M{R) space that can be 
reached by hybrid quark-nuclear stars are the same as those that can be reached by 
hadronic matter stars, once moderately exotic phenomena such as kaon condensation 
or hyperon production are allowed Observation of an object with a small radius 
(i? ~ 7 to 10 km at M 1.4 Mq) would rule out simple nucleonic matter, but would 
not favor quark matter over the exotic forms of hadronic matter. There are regions 
of parameter space (very small M and R, for example) that can only be reached 
by pure quark matter objects, which only exist if quark matter is absolutely stable 
(pc = 0). Other regions of the parameter space are inconsistent with both hadronic 
and quark matter. 

Demonstrating the presence of color- superconducting quark matter via M{R) 
measurements appears to us to be very difficult. However, this is not for the naive 
reason that "the color superconductivity contribution to the pressure is suppressed 
by C(A^//i^)". As previously noted IHJO], the leading fi'^ contribution is mostly 
canceled by the bag constant, so the subleading A^/i^ term is potentially important. 
However, in practice it is not detectable. Firstly, the bag model equation of state 
depends on the color superconducting gap A via the linear combination — 4A^, 
so an accurate determination of nis at high density would be needed to expose 
the presence of a nonzero A. Secondly, the dependence of the M{R) relation on 
ml — 4A^ is not particularly strong (see Fig. [7j) when a physical parameter pc is 
assumed known, rather than the bag constant. 

If color-superconducting quark matter is to be found in compact stars, it 
seems more likely that it will be detected via its effects on transport properties. 
Color superconductivity drastically alters these, and possible signatures are being 
actively investigated. These include cooling (sensitive to heat capacity and neutrino 
emissivity and opacity) jl^j ? T-mode spin-down (sensitive to bulk and shear viscosity) 
jlH], and glitches (sensitive to superfluidity and rigid structures) [T3] . 
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